Theory of carrier mediated ferromagnetism in dilute magnetic oxides 
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We analyze the origin of ferromagnetism as a result of carrier mediation in diluted magnetic oxide 
semiconductors in the light of the experimental evidence reported in the literature. We propose 
that a combination of percolation of magnetic polarons at lower temperature and Ruderman-Kittel- 
Kasuya-Yosida ferromagnetism at higher temperature may be the reason for the very high critical 
temperatures measured (up to ~ 700K). 

PACS numbers: 75.50.Pp, 75.30.Hx, 75.10.-b 
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I. INTRODUCTION 

Semiconductors doped with magnetic ions are being 
studied in an effort to develop spintronics, the new kind 
of electronics that seeks to exploit, in addition to the 
charge degree of freedom as in the usual electronics, also 
the spin of the carriers [ij. The first so-called diluted 
magnetic semiconductors (DMS) were II- VI semiconduc- 
tor alloys like Zni—a, Mna;Te and Cdi-a^Mn^^Te [i| orig- 
inally studied in the 1980s. These materials are either 
spin glasses or have very low ferromagnetic (FM) critical 
temperatures Tq (~ few K) [sj and are, therefore, in- 
adequate for technological applications which would re- 
quire FM order at room temperature. More recently, the 
Mn doped III-V semiconductors Ini^^jMn^jAs [1, Q and 
Gai_a;Mn2,As ^ {7^ showed ferromagnetism at a much 
higher temperature, thanks to the development of molec- 
ular beam epitaxy (MBE)-growth techniques. The cur- 
rent high Tc record of 173K achieved in Mn-doped GaAs 
by using low temperature annealing techniques [1, 
is promising, but still too low for actual applications. In 
all these materials, ferromagnetism has been proven to 
be carrier mediated, a necessary property for spintronics 
since this enables the modification of magnetic behav- 
ior through charge manipulation. This has motivated 
a search for alternative spintronics materials with even 
higher Tc and carrier mediated FM. In this direction, 
dilute magnetic oxides such as magnetically-doped 
TiOa [il, ZnO [H, and SnOz [ll, could represent this 
alternative with reported Tqs above room temperature 
and as high as - 700K [H. 

The general formula for these oxide based dilute mag- 
netic semiconductors (0-DMS) is 



(1) 



where A is a non-magnetic cation, M is a magnetic 
cation, S is the concentration of oxygen vacancies which 
depends on the growth conditions, and n = 1 or 2. Car- 
riers, usually electrons, are provided by the oxygen va- 
cancies that are believed to act as shallow donors [l^, 
[iTj . This is in contrast to III-V semiconductors, like 
Gai_a;Mn3;As, where the carriers (holes in this case) are 
provided by the magnetic impurities themselves which 
act also as donors (or acceptors). More oxygen vacan- 



cies and, thus, niore carriers, are produced at low oxygen 
pressure [H, [l^. In the process of doping with mag- 
netic ions, usually of different valence than the ion they 
substitute for, oxygen or A vacancies are also introduced 
to maintain charge neutrality. However, these oxygen 
vacancies have been found not to contribute directly to 
the electrical conductivity of the system On the 

contrary, the resistivity increases by orders of magnitude 
upon doping [13, [l^ . 

There is currently no consensus on the origin of fer- 
romagnetism in 0-DMS, in particular, whether it is an 
extrinsic effect due to direct interaction between the local 
moments in magnetic impurity clusters (or nanoclusters) 
or is indeed an intrinsic property caused by exchange 
coupling between the spin of the carriers and the local 
magnetic moments. This is a very important issue be- 
cause spintronics requires the carriers to be polarized and 
this would only be guaranteed if ferromagnetism is in- 
trinsic. Experimental evidence for carrier-mediated FM 
in 0-DMS is not yet conclusive. In the much studied 
Co-doped Ti02 in the rutile phase, anomalous Hall effect 
(AHE) 18] and electric field induced modulation of mag- 
netization (by as much as 13.5%) ^l| have been observed, 
arguing for carrier-mediated FM. However, AHE was also 
measured in a sample with magnetic clusters, casting 
doubts about the conclusions that can be drawn from 
AHE data [13]. Nevertheless, room temperature FM has 
also been reported in cluster free films grown or annealed 
at high temperatures ~ 900°C IB'] by pulsed laser depo- 
sition, and in nanocrystalline films grown in conditions 
that preclude the formation of metallic cobalt [1^. An- 
other piece of evidence in favor of carrier-mediated ferro- 
magnetism is the observation that the magnetic field de- 
pendence of ferromagnetic circular dichroism is in good 
agreement with those of magnetization and anomalous 
Hall effect [IJ] . More recently. X-ray photoemission spec- 
troscopy measurements have suggested strong hybridiza- 
tion between carriers in the Ti 3d band and the localized 
t2g states of Co2+ [H]. 

Many reports have raised serious doubts on the mag- 
netism of magnetically-doped ZnO [2^ since the results 
are very sensitive to sample preparation. However, it 
has been pointed out [1^, |23] that the lack of ferromag- 
netism in some samples can be the result of too low a 
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density of carriers. Indeed, films of Zno.75Coo.25O pre- 
pared in low oxygen partial pressure (< 10~^ Torr), a 
condition that should increase the density of electrons 
from oxygen vacancies, were found to be ferromagnetic 
at room temperature [28| . These samples exhibited per- 
pendicular magnetic anisotropy [H, Hsj and no segre- 
tion effects were significant. In the same direction, dop- 
ing Zno.98Coo.02O with small amounts of Cu also en- 
hanced the system's ferromagnetism [s^]. The system- 
atic variation of magnetism in doped ZnO as a function 
of the magnetic dopant has been explored in Ref. [l9[, 
where room-temperature ferromagnetism has been found 
in films doped with Sc, Ti, V, Fe, Co, or Ni but not with 
Cr, Mn, or Cu. The same group also reported an increase 
in the magnetic moment per Co ion as a function of the 
reduction of oxygen pressure (equivalent to an increase 
of the density of carriers) [l^ . Optical magnetic circular 
dichroism, one of the tests proposed as a signature of di- 
luted ferromagnetism [2]| , has also been measured at low 
temperatures |3l| . and at room temperature {s^]. Very 
recently, by analyzing the controlled introduction and re- 
moval of interstitial Zn, carrier-mediated ferromagnetism 
in Co-doped ZnO has been demonstrated [33l |. 

The case of magnetically-doped Sn02 (with Tc ~ 
600K seems to be different from the previous two 

materials in that the parent compound is highly conduct- 
ing [sl] (though still transparent) and that some doped 
samples have shown an extremely lar ge m agnetic moment 
(> 7/iB) per magnetic impurity l33|. 

It thus appears that different doped magnetic oxides 
may very well have different underlying mechanisms lead- 
ing to the observed ferromagnetism - in particular, there 
may very well be several competing FM mechanisms in 
play in 0-DMS, a novel theme we develop further here. 

In this article, we address the origin of ferromagnetism 
as mediated by carriers in 0-DMS analyzing the exist- 
ing experimental evidence for magnetic and transport 
properties, in particular, in Co-doped Ti02. Our goal 
is to develop a theory for 0-DMS ferromagnetism (as- 
suming it to be mediated by band or impurity-band 
carriers) in analogy with the better understood III-V 
Mn-doped DMS materials 0, Hi, [13, S [3i|. Magnet- 
ically doped III-V semiconductors are ferromagnetic for 
a wide range of carrier concentrations, from the insulat- 
ing to highly conducting regimes, where different mech- 
anisms are expected to prevail. It has been proposed 
that for high enough values of the local exchange be- 
tween the carriers and the magnetic ions, an impurity 
band is formed [13, |4l|. In a highly insulating system, 
the Fermi level is well below the mobility edge of the im- 
purity band. In this regime, ferromagnetism has been 
explained as the result of percolation of bound mag- 
netic polarons [4^ . This mechanism is consistent with 
the concave shaped magnetization versus temperature 
M{T) curves observed, for instance, in (In,Mn)As Q. 
In the more conducting samples, itinerant carriers would 
mediate ferromagnetism via a Ruderman-Kittel-Kasuya- 
Yosida (RKKY) mechanism whose sign fluctuations, and 



the associated frustration effects, are suppressed due to 
the low density of carriers [43| . arising possibly from 
heavy compensation and/or existence of defects. 

Given that the carriers in magnetic oxides reside in 
an insulating impurity band, there are essentially three 
kinds of carrier-mediated magnetic exchange interac- 
tions which could potentially lead to intrinsic carrier- 
mediated ferromagnetism: double exchange (similar to 
the situation in manganites) (4^ in the impurity band, 
bound magnetic polarons percolation (similar to the situ- 
ation in insulating diluted magnetic semiconductors, eg. 
Gei_j;Mnj;, Ini-ajMua-As, and low- a; Ga.i-xMi^x^'^) j43 |. 
and indirect (RKKY) exchange coupling mediated by free 
carriers (similar to the situation in the optimally doped 
X « 0.05 high-Tc metalhc Gai_^Mn^As) glj. The dou- 
ble exchange mechanism gives, at the low carrier density 
of magnetic oxides, Tc proportional to the carrier den- 
sity; therefore, critical temperatures exceeding room tem- 
perature (~ iQQK) are essentially impossible within this 
model. In general, the carrier-mediated indirect RKKY 
exchange mechanism applies only to free carriers in an 
itinerant band, and therefore may be thought to be ruled 
out for doped magnetic oxides which are insulating im- 
purity band materials. But we propose here that the 
'standard' RKKY mechanism may very well be play- 
ing a role in the magnetic oxides, particularly at high 
temperatures, where a large number of localized impu- 
rity band carriers will be thermally activated (due to the 
small values of the carrier binding energies in the mag- 
netic oxides) either within the impurity band or to the 
conduction band, effectively becoming itinerant free car- 
riers (albeit thermally activated ones) which can readily 
participate in the indirect RKKY exchange between the 
localized impurity magnetic moments. If the exchange 
coupling between local moments and the thermally acti- 
vated carriers is sufficiently large, then this mechanism 
could explain the very high Tc observed in magnetic 
oxides. We emphasize that Tc also depends on other 
parameters including the activated carrier density and 
the magnetic moment density. At low temperature, the 
RKKY mechanism must freeze out since thermal acti- 
vation from the impurity band is no longer operational 
and one must therefore have a complementary mecha- 
nism to provide carrier-mediated ferromagnetism. We 
believe that the bound magnetic polaron (BMP) mecha- 
nism is the plausible low temperature magnetic ordering 
mechanism, but it cannot explain the high (room tem- 
perature or above) claimed Tc of 0-DMS unless one uses 
completely unphysical parameters. 

Our new idea in this work is, therefore, the suggestion 
that the intrinsic carrier-mediated ferromagnetism lead- 
ing to high Curie temperatures is plausible in doped mag- 
netic oxides only if two complementary magnetic mech- 
anisms (i.e. the bound magnetic polaron percolation 
at low temperatures and the indirect RKKY exchange 
mechanism in the presence of substantial thermal acti- 
vation of carriers in the impurity band) are operating in 
parallel. We find that no single carrier-mediated mech- 
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anism, by itself, can account for the observed high Tc 
in doped magnetic oxides. These considerations apply 
mainly to Ti02 as high-T resistivity measurements indi- 
cate existence of substantial thermally activated carrier 
population [l^[ll,|4^. Of course, the possibility that the 
observed ferromagnetism in magnetic oxides arises from 
completely unknown extrinsic mechanisms (e.g. cluster- 
ing near the surface) or from non-carrier-mediated mech- 
anisms cannot be ruled out at this stage. Only more 
experimental work, perhaps motivated by our theoreti- 
cal considerations presented herein, can provide defini- 
tive evidence in favor of one or more FM mechanisms in 
0-DMS. We mention in this context that even for much 
better-understood III-V DMS systems there are still de- 
bates in the literature regarding the precise role of the 
RKKY mechanism [46j . 

This paper is organized as follows: in Sec.|TTl the bound 
magnetic polaron and the RKKY model are introduced. 
In Sec. Illli magnetic and transport properties of dilute 
magnetic oxide semiconductors, mainly Co-doped Ti02, 
are summarized and analyzed in the light of our pro- 
posed combined model. Sec. IIVI presents a discussion of 
our model and some of the alternatives suggested in the 
literature. We conclude in Sec. |Vl 

II. MODEL 

The general Hamiltonian that describes the physics 
of diluted magnetic semiconductors is (see, for instance, 
Ref. [13) 

+ J d^x^ JS, 'ii{x)ac,,p^p{x) 4S^{x-R^) 

ia/3 

+ J2Jd{R^-RJ)SiSJ (2) 

i,j 

where m is the relevant effective mass, Vl is the peri- 
odic lattice potential, Vr is a potential arising from dis- 
order (magnetic and non-magnetic) in the lattice, W is 
a Coulombic potential arising from the oxygen vacan- 
cies that act as shallow donors, J is the local exchange 
(Hund's like) between the carrier spin and the magnetic 
impurities moments, and Jd is a direct exchange between 
the magnetic im pur ities spins (which is ferromagnetic in 
Co-doped Ti02 [4g|). is the impurity spin located at 
R-i, fct./S represents the Pauli matrices with spin indices 
a and (3, and Qq is the unit cell volume. We are neglect- 
ing the electron-electron interaction which we expect to 
be very small due to low carrier density in the diluted 
magnetic oxides. In the following, we assume the direct 
exchange term to be effectively included in the local ex- 
change term. Only at very low carrier densities and in 



the case of antiferromagnetic Jd, which is not the case of 
Co-doped Ti02, the direct exchange could compete with 
the carrier- mediated FM and cause frustration [i^ . 

The Hamiltonian in Eq. ^ is extremely complex to 
solve exactly and, consequently, we integrate out the elec- 
tronic degrees of freedom and simplify the problem by 
considering only the term in the Hamiltonian that dic- 
tates the interaction between the spin of the carriers and 
the magnetic moments 

Hi^J2 J"-o s(R,) . (3) 

i 

We, therefore, model the problem with a minimum set 
of parameters that effectively include other interactions 
in the system. This simplified Hamiltonian is solved in 
two complementary cases: (i) localized carriers (bound 
to oxygen vacancies by the interaction W), and (ii) itin- 
erant carriers. When the carriers are localized, they form 
bound magnetic polarons that percolate at Tc (43 . |50| . 
When carriers are delocalized in the conduction (or va- 
lence) band, they can mediate ferromagnetism through 
the RKKY mechanism. In the following, we introduce 
these two approaches, which are capable of producing 
carrier-mediated FM provided the carrier density is much 
lower than the magnetic impurity density. 

A. Percolation of bound magnetic polarons 

Bound magnetic polarons (BMP) are the result of 
the combination of Coulomb and magnetic exchange in- 
teractions [5l|. Carriers are localized due to electro- 
static interaction with some defect (i.e. the magnetic 
impurity in III-V semiconductors, and the oxygen vaca- 
cies in the magnetic oxides) with a confinement radius 
flB = e(m/m*)a, with e the static dielectric constant, 
TO* the effective mass of the polaron, and a — 0.52 A 
the Bohr radius, ob is larger (smaller) for shallower 
(deeper) defect levels. We assume the bound electron (or 
hole) wave-function has the isotropic hydrogen-like form 
vl'(r) ~ — i= exp(— r/ae) so that the exchange [given in 

Eq. ([3])] between the magnetic impurities and the carrier 
decays exponentially with the distance r between them 
as ^ exp(— 2r/aB). At a certain T, the radius of the 
polaron Rp is given by the condition 

ksT = \J\ (ao/as)' Ss exp(-2i?p/aB) , (4) 
where /cb is the Boltzmann constant, leading to 

i?p(T) = (aB/2) In J| {ao/ajif /fesT) , (5) 

where we can see that at the high temperature disor- 
dered phase there are no polarons, which start forming 
at ksT ~ sS\J\ (oo/aB)'^, and then grow as temperature 
is lowered. The impurity spins that are a distance r < Rp 
from a localized carrier tend to align with the localized 



4 



carrier spin. The polarized magnetic impurities form, in 
turn, a trapping potential for the carrier such that a finite 
energy is required to flip its spin. The magnetic polarons 
are well defined non overlapping isolated entities only at 
low carrier densities and sufficiently large temperatures. 
The size of the polarons increases as temperature de- 
creases, eventually overlapping with neighboring BMPs. 
This overlap causes the alignment of their spins, there- 
fore forming FM clusters. The FM transition takes place 
when an 'infinite cluster' (of the size of the system) is 
formed, i.e. when the BMP percolation occurs. This sce- 
nario has been studied in Ref. [i^ in the context of III-V 
semiconductors. There it is shown that the bound mag- 
netic polaron model can be mapped onto the problem of 
percolation of overlapping spheres. The model proposed 
is valid in the low carrier density regime nca% << 1 and 
when the density of magnetic impurities Ui is larger than 
the density of carriers ric- Under these conditions, each 
carrier couples to a large number of magnetic impurities, 
as shown in Fig. [TJ 

In order to calculate the BMP percolation critical tem- 
perature Tq^'^'^ we first need to estimate the maximum 
temperature at which two magnetic polarons a distance 
r from each other are still strongly correlated T2p(r). By 
estimating the number of impurities which interact with 
both polarons ^45| , this temperature is found to be given 

by 



kBT2p{r) ^ aBy/rn^sSJ 



exp(-r/aB). (6) 



As the temperature is lowered more and more polarons 
overlap with each other until a cluster of the size of the 
sample appears. The critical polaron radius at which 
this happens can be calculated as the percolation radius 
in the problem of randomly placed overlapping spheres. 
This problem has been solved numerically [54I giving 



0.86/^ 



(7) 



Substituting this distance in Eq. [S] gives the FM transi- 
tion temperature 



B-Lc 



Q T I \ I 3 ^l/3 Ui 
sSJ — (annc) \ — exp 
\aB/ V 



0.86 



.(8) 

The magnetization is due to the magnetic ions (since 
rii > Tic) in the percolating cluster. It is given by 



(9) 



where V{y) is the infinite cluster's volume in the model 
of the overlapping spheres and 



0.86 + (a|ne)i/3ln^ 



(10) 



is the maximum distance between correlated polarons at 
T. Therefore, the shape of the magnetization curve only 




FIG. 1: Schematic representation of magnetic percolation 
in oxide based dilute magnetic semiconductors. The solid 
squares represent the oxygen vacancies where an electron, rep- 
resented by an arrow, is localized. The gray circles represent 
the extension of the electron wave- function. Magnetic impu- 
rity spins are represented by small arrows whose orientation is 
established by antiferromagnetic exchange coupling with the 
localized carrier. 



depends on 



The volume of the infinite cluster 



is strongly suppressed as temperature is increased pro- 
ducing a concave shaped M{T) curve as opposed 
to the Weiss mean-field like convex Brillouin function 
shape expected for highly conducting samples. The con- 
cave shape is less pronounced for larger a^ric [l^l, that 
would eventually lead to a change of shape in the itin- 
erant carrier regime, beyond the limit of applicability of 
this theory. Interestingly, however, the BMP percola- 
tion theory developed in the strongly localized insulating 
regime smoothly extrapolates to the mean-field RKKY- 
Zener theory (discussed below in SecHTb) in the itinerant 
free-carrier metallic regime [4^ l47j , giving us some addi- 
tional confidence in the potential validity of the compos- 
ite percolation-RKKY model we are proposing here for 
ferromagnetism in the magnetic oxides materials. 

Conduction in a bound magnetic polaron system oc- 
curs via an activated process. Thus, p exp{AE / kBT) 
where AE is the activation (or binding) energy. AE de- 
pends on the binding energy of the carrier due mainly 
to electrostatic interaction with donors or acceptors, and 
in a minor degree, on the polarization due to magnetic 
exchange Eq. ([3]). The effect of an applied magnetic field 
in this system is to align non-connected polarons in such 
a way that the polarization part of the binding energy 
gets suppressed. As a result, the magnetoresistance is 
negative and depends exponentially on the field [53,]. It 
should be emphasized that the BMP picture and the as- 
sociated magnetic percolation transition is only valid in 
the strongly localized carrier regime where each BMP 
(i.e. each carrier) is immobile, and a percolation transi- 
tion in the random disordered configuration makes sense. 
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B. RKKY model for itinerant carriers 

As already mentioned, carriers in magnetic oxides are 
donated by oxygen vacancies that act as shallow donors. 
The carrier binding energy has been estimated for differ- 
ent samples of undoped Ti02 to be '--^ 4 meV [16^ , ~ 14 
meV [13 and ^41 meV jl^l- These are smaller than the 
Mn acceptor binding energy of GaAs, which is of the or- 
der of 0.1 eV dl]. It is then possible that a large fraction 
of the electrons are promoted (i.e. thermally activated) 
above the mobility edge in the impurity band at temper- 
atures lower than T^°'^'^ and that the localization picture 
just described cannot explain ferromagnetism by itself. 
This thermal activation of carriers will be extremely ef- 
fective at higher temperatures, particularly since k^Tc in 
this system could be substantially higher than the elec- 
tron binding energy or, equivalently, the activation en- 
ergy. The assumption of strict BMP localization fails in 
this situation as the carriers get activated into highly mo- 
bile states at higher temperatures, as manifested by the 
experimentally observed enhancement in the conductiv- 
ity at higher temperatures. The BMP formation will be 
suppressed in this situation since the thermally activated 
band carriers are effectively free or itinerant at elevated 
temperatures. 

Itinerant carriers coupled to local moments by Eq. ([3]) 
lead to the well-known Zener-RKKY mechanism of indi- 
rect magnetic interaction between the magnetic impuri- 
ties. This model gives an effective exchange Jpfp between 
magnetic moments in the lattice of the form [SQ . Is^ 



Joff oc 



2kpRij cos(2/ci?i?y ) — sm(2kpRij) 



(11) 



where Rij — Ri — Rj (i.e. the distance between magnetic 
impurities) and kp is the Fermi momentum. Jeff oscil- 
lates in sign with distance leading in general to compli- 
cated magnetic configurations [58i] , frustration and glassy 
behavior. However, in the limit of very low density of 
carriers relevant here {uc << Ui), kpr and Jeff < 0, 
namely, the RKKY interaction is always ferromagnetic. 
The RKKY ferromagnetism phase diagram has recently 
been calculated [4^ . 

We study the RKKY interaction in the limit of non- 
degenerate carriers within mean-field theory (therefore, 
all disorder in the lattice is neglected). The interaction 
between carrier spins and local moments can be then de- 
scribed as a self-consistent process in which carrier spins 
see an effective magnetic field produced by the local mo- 
ments (which are considered classical) 



-°cff — 

gc^J-B 



(12) 



and the local moments see an effective field produced by 
the carrier spins 



B 



(13) 



The response of the impurity spin, which follows Boltz- 
mann statistics, to B^g is [59|] 



(5,) = SBs 



knT 



= SBf 



V kBT 



where 



^ , , 2s + 1 , 2s + 1 1 1 



(14) 



(15) 



is the Brillouin function. In the nondegenerate limit, 
the carrier spin distribution is not affected by the Pauli 
exclusion principle and, therefore, it is determined also 
by Boltzmann statistics rendering 



knT 



sJalni{Sz 
kT,T 



(16) 



Combining Eqs. ([Tl)) and we get a self-consistent 

equation for the impurity spin 



SBs 



Jap 
knT 



sBs 



Jalni{Sz 
keT 



(17) 



When the effective magnetic fields are small, as it is close 
to the magnetic transition, we can perform the expansion 
of Bs{y) for small values of y 



Bsiv) 



1 



3s 



J/ + 0(y3), 



(18) 



which, applied to Eq. I17[ gives the critical tempera- 
ture lei 



- Jail ^{S+l){s + l). (19) 



The magnetization is mainly due to the local magnetic 
moments {Sz) (because rii > Uc) and is given by the 
self-consistent Eq. [171 The resultant dependence with 
temperature is concave for low values of ric/ni and convex 
mean-field-like for Tic/rii > 0.2 (47| . 

For higher carrier density, the carriers form a degener- 
ate gas and the procedure followed in Ref. [1^ would be 
more appropriate. In this reference, the RKKY model is 
solved taking into account the disorder in the lattice. The 
conductivity is included via the mean free path (MFP), 
which produces a cutoff in the RKKY coupling reach. Tc 
has been found to depend very strongly on the relation 
between ric, Ui and the MFP. MFP can be made larger 
by improving sample quality. In contrast with mean field 
treatments of the same model (which neglect disorder) 

that predict Tc oc n]/^ [6l|, Tc is enhanced and later 
suppressed by increasing ric- This Tc suppression is due 
to the sign oscillations of the RKKY interaction at larger 
carrier densities [see Eq. (fTTj) ] that lead to magnetic frus- 
tration and spin glassy behavior. It is then proposed [6^ 
that Tc improvement can arise both from increasing ric 
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and the MFP. The magnetization curves M(T) are con- 
vex in the highly conducting Hmit (large MFP) and con- 
cave for more insulating systems (small MFP) [i^, |43| 
appropriate for dilute magnetic oxides. This is quali- 
tatively the same result as given by the non-degenerate 
approach and the BMP model. In this way, the localized 
and itinerant carrier models predict the same behavior 
for M{T) in both the high and low carrier density limits, 
leading to correct qualitative predictions even beyond the 
applicability range of the two models. 

Calculations of dc-resistivity in diluted systems within 
DMFT approximation have rendered a negative magne- 
toresistance MR, that peaks at Tc but is appreciably 
smaller than the MR of magnetic ordered lattices such 
as manganites (63| . 

C. Double exchange model 

In the double exchange model [l^l , a large J forces the 
spin of the carrier to be parallel to the local magnetic 
moment in such a way that the kinetic energy of the car- 
riers hopping between magnetic sites is minimized when 
the magnetic ions are ferromagnetically ordered. This 
model was proposed for manganites, where ferromag- 
netism and metallicity usually come together. The fer- 
romagnetic critical temperature in the double exchange 
model is proportional to the density of carriers and the 
bandwidth |64| , with very small values in the low density 
limit appropriate to 0-DMS. Therefore, double exchange 
is not a suitable mechanism to explain the high critical 
temperature of dilute magnetic oxides, and is hence not 
considered further in this work. However, in the strong 
coupling regime, where the effective exchange coupling J 
is large {J >> t, Ep where t is the carrier bandwidth), 
a double exchange mechanism may be appropriate for 
low-Tc DMS materials as discussed in Refs. [io, [6^ . 

III. COMPARISON TO EXPERIMENT 

We argue here that, at high enough temperatures com- 
pared to the carrier binding energy, carriers are thermally 
excited to conduction (or valence) bands from the impu- 
rity band becoming itinerant and mediating ferromag- 
netism by an effective RKKY mechanism. The carrier 
binding energy AE can be estimated by fitting the re- 
sistivity curves to p = poex.p{AE/kBT). The result is 
unfortunately sample and composition dependent. For 
the undoped Ti02, estimations go from ^ 4 meV 16] to 

41 meV dll. In Ref. [l3|, as in Ti02 is estimated 
to be 15 A from the observation of an insulator to metal 
transition upon doping at ric ^ 5 x 10^^ cm~'^. This 
corresponds to AE = 14 meV using the static dielectric 
constant e = 31 [l3| and m* — m. We show these results 
in Fig. [21 together with the activation energy values cor- 
responding to Co-doped Ti02 films [53|. When doped 
with magnetic ions, the resistivity increases by orders of 




%Co 

FIG. 2: Activation energies (closed symbols) and estimated 
localization radius ob (open symbols) for undoped and Co- 
doped Ti02. For these estimations, the static dielectric con- 
stant e = 31 has been used. Circles are results from Ref. [53 ]. 
squares from Ref. ^V?\, and diamonds from Ref. p^ . 



magnitude proportionally to the density of impurities, 
its slope is dramatically enhanced, and does not show 
the insulator to metal transition [ij, [l^ [11] observed in 
the undoped samples. The overall increase of the resis- 
tivity is due to the scattering from charged impurities, 
as they produce strong disorder in the system. However, 
the resistivity due to the impurity scattering is not tem- 
perature dependent on the exponential scale of carrier 
activation. Therefore, AE can be calculated the same 
way as for undoped samples. The results for a particular 
series of Co-doped Ti02 ^54] show that the activation en- 
ergy increases upon doping (see Fig. (2). We expect these 
values to vary widely from sample to sample, as observed 
for the undoped compound. (This is also consistent with 
the current experimental situation where the observed 
Tc in various nominally similar Co-Ti02 samples varies 
widely from sample to sample.) 

As the density of free carriers depends exponentially 
on the impurity energy level (or the activation energy), 
a significant part of the bound impurity band electrons 
can be promoted above the mobility edge at tempera- 
tures lower than Tc- This is particularly true in doped 
magnetic oxides, where the claimed Tc (~ 400 — lOOOK) 
is so high that the /cbTc >> AE regime is easily reached 
producing a high density of thermally activated mobile 
band carriers at fceT > AE. 

The estimation of the critical temperature depends 
on many variables: (i) the density of carriers Uc, (ii) 
the density of magnetic ions n^, (iii) the exchange cou- 
pling J and, (iv) in the bound magnetic polaron theory, 
the localization radius ae. These parameters are either 
not precisely known or present a huge sample-to-sample 
variation: (i) The density of carriers ric has been mea- 
sured by Hall effect giving results ranging from 10^* to 
10^^ cm"'^ [4^, as the density of oxygen vacancies varies 
greatly with the growth conditions, (ii) The absolute 
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density of magnetic ions is known (a; = 0.1 corresponds 
to ~ 3 X 10^^ cm~'^ in Ti02) but we cannot estimate 
how many of those are magnetically active and, therefore, 
relevant to long-range carrier-mediated ferromagnetism. 
In fact, the magnetic moment per magnetic ion has a 
strong dependence on sample characteristics [50]. This 
could be due to different effective values of magnetically 
active impurities. In particular, for Co-doped Ti02, 5* is 
usually close to the low spin state S = 1/2. (iii) J is not 
known in general though there are some estimates which 
place its value above 1 eV for ZnO [H^l- (iv) ob calcu- 
lated from the activation energies are shown in Fig. [5] 
(right axes). Note the strong variation of ob from sam- 
ple to sample in the undoped case. This dramatically af- 
fects the estimates of the critical temperature within the 
polaron percolation model as shown in Fig. [31 as Tq'^^ 
depends exponentially on a^ric- 

Estimations of Tc from Eqs. ([5]) and depend 
strongly on all these unknown parameters and, conse- 
quently, Tc can vary from tens to hundreds of Kelvin 
upon tuning their values. This is illustrated in Fig. [31 
The local moment is taken to be 5' = 1, however, the 
moment per magnetic ion is another quantity that is very 
sample dependent. A value of = 3 x 10^^ cm~^ has 
been used. TjS"'^'^ [Eq. ([5])] increases with the magnetic 
impurity density as ^/ni while JJ^J^kky jgq^ (O] only 
depends on the product y/niJic- Note that both esti- 
mates of Tc are the same order of magnitude and close 
to experimental data. However, due to the strong depen- 
dence of Tc on unknown parameters, we should not use 
the calculated value of Tc as the sole criterion to elu- 
cidate the applicability of a particular model. Rather, 
we should look at other evidence given by experiment 
such as trends in magnetization curves and magnetore- 
sistance. We note, however, that we cannot rule out the 
possibility that the strong variation in Tc between differ- 
ent experimental groups (or even from sample to sample 
in the same group) arises precisely from the variation of 
the sample parameters n^, ric, ob, etc. which will in- 
deed lead to large Tc variation! Further experiments are 
clearly needed to settle this important issue. 

The measured magnetization versus temperature 
curves M[T) usually present a very constant signature 
within a wide range at low temperatures (T^ . [3, lisl [6^ . 
though in some cases, a concave shape is observed [2^, 
l68j . As explained in Sec. [Til tbis concave shape is ex- 
pected for insulating systems within the polaron percola- 
tion approach and for low carrier density systems within 
the RKKY approach. 

Magnetoresistance is very sample and composition de- 
pendent and there are not many reports in the literature. 
In general, it is only significant at temperatures much 
lower than Tc [ll,[l4, 69]. For doped ZnO, MR shows dif- 
ferent signs for different dopants and temperatures [69| . 
though these results may not be relevant as these samples 
were not ferromagnetic. Co-Sn02 displays positive MR 
that becomes negative with increasing T [IJ]. Anatase 
Co-Ti02 has positive MR for T < 5 K (60% at T = 2K 
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FIG. 3: (a) Estimation of ferromagnetic Tc for the theory 
of bound magnetic polarons Eq. (|8]) and (b) for the RKKY 
model for itinerant carriers Eq. (I19p . Tc on the right y-scale 
has been calculated using J =1 eV. S = 1, ao = 3.23 A 
and Tii = 3 X 10^^ cm"'^ (equivalent to a; = 0.1 for doped 
TiOa). 

and H — 8T) [12], while rutile Co-Ti02 shows negative 
MR up to almost room temperature (maximum ~ —0.4% 
at r = lOOK and = 8T) [11]. The latter also presents 
an increase of MR for the least resistive samples, the ones 
with higher carrier density. 

As explained above, the polaron percolation model 
can explain a negative MR, whose value depends on the 
change of {Sz) with the applied magnetic field. At low 
temperatures, {Sz) should already be saturated, there- 
fore, MR is bound to be small, consistent with results 
in 0-DMS. On the other hand, both the percolation 
and RKKY models imply a maximum MR at Tc when 
spin fluctuations and susceptibility are very large and a 
small magnetic field can affect magnetization dramati- 
cally. However, this behavior has not been found in O- 
DMS. 



IV. DISCUSSION 

In this work we are trying to understand the physi- 
cal origin of ferromagnetism in magnetically doped ox- 
ide materials (e.g. Co-Ti02) by assuming the underly- 
ing magnetic mechanism to be carrier-mediated, qualita- 
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tively similar to that operational in the widely studied 
DMS systems such as Gai_a;Mna;As. We use an effec- 
tive model of a local exchange coupling between the lo- 
calized magnetic moments (i.e. Co) and the localized 
carriers (since the system, e.g. Co-Ti02, is an insulator 
at T — 0). We first show that the very high Tc (well 
above the room temperature) reported in the literature 
cannot be understood entirely within the BMP perco- 
lation picture by itself without fine-tuning the effective 
parameters (e.g. exchange coupling, carrier and mag- 
netic moment densities) to an unreasonable degree. We 
therefore disagree with the recent attempts [50] in the lit- 
erature to attribute FM in 0-DMS entirely to BMP per- 
colation. In general, the BMP percolation theory, which 
is very natural for the insulating magnetic oxides and 
has already been invoked in the literature [l^ [H, [soj . 
leads to low values of Tc for any reasonable assump- 
tions about the system parameters. The observed convex 
shape of the magnetization curves also argues against an 
entirely BMP percolation model for magnetic oxide ferro- 
magnetism. In addition, for Ti02, the high temperature 
resistivity measurements indicate the existence of a sub- 
stantial thermally activated carrier population consistent 
with RKKY. 

We therefore suggest a composite scenario by taking 
into account the relatively small carrier binding energy 
AE in these materials. We propose that the low temper- 
ature BMPs in the system (in the ferromagnetic state) 
become mobile (instead of going into a disordered para- 
magnetic state through the percolation transition), and 
these mobile carriers then produce long-range ferromag- 
netic coupling among the magnetic impurities through 
the standard RKKY-Zener mechanism, which could ex- 
plain the experimentally observed high-temperature fer- 
romagnetism. Wide variations in effective carrier (n^) 
and magnetic moment (n^) densities due to different 
growth conditions could easily accomodate the observed 
large variation in Tq among different groups. 

We emphasize that our model of combined polaron per- 
colation (at low temperatures, where the carriers are still 
bound) and RKKY-Zener (at high temperatures, where 
the carriers are activated into mobile states) mechanism 
seems to be the only reasonable possibility for the fer- 
romagnetism with high Tc of the doped magnetic ox- 
ides assuming that the ferromagnetism is both intrin- 
sic and carrier-mediated. This combined model applies 
mainly to TiOo where high temperature resistivity is ac- 
tivated [l^, [H, H^. On the other hand, recent mea- 
surements on Co-doped ZnO seem to rule out RKKY in 
favor of an exclusive BMP mechanism 33]. There are 
obviously other possibilities for the origin of ferromag- 
netism, which are less interesting, but not necessarily 
less plausible. First, the oxide ferromagnetism (in Co- 
Ti02, Co-ZnO, etc.) could be entirely extrinsic (and, 
some of it most likely is) , arising from contaminants (as 
happened, for example, in the ferromagnetism of cal- 
cium hexaborides) or from magnetic nanoclustering . 
The Co nanoclusters in Co-Ti02, for example, could pro- 



duce long-range ferromagnetic order through dipolar cou- 
pling. A second intrinsic possibility for oxide ferromag- 
netism, increasingly gaining ground on the basis of re- 
cent first principles band structure calculations, is that 
the ordering in insulating samples or regions arises from 
the standard superexchange interaction between the lo- 
calized magnetic impurities (e.g. Co in Co-Ti02) ^48|. 
These local interactions would still coexist with long- 
range carrier-mediated ferromagnetism. Other calcula- 
tions (70j find a, p — d hopping mechanism, not related to 
free carriers but to clustering of the magnetic ions. These 
possibilities certainly cannot be ruled out, particularly 
in ordered magnetic oxides where such band structure 
supercell calculations would, in principle, be applicable. 
But the first principles calculations completely neglect 
disorder, and cannot really explain the wide variation in 
magnetic properties (e.g. Tc) seen in different labora- 
tories. In addition, none of these posibilitics (superex- 
change, p-d hopping mechanisms, or extrinsic nanoclus- 
ter ferromagnetism) are carrier-mediated in the sense of 
DMS materials, and as such are out of scope for our in- 
terest. 

What we have discussed in this work is that, within 
the assumption of an intrinsic carrier-mediated ferro- 
magnetic mechanism (related to the mechanisms oper- 
ational in DMS materials), the most likely origin for fer- 
romagnetism in DMS oxides is a combination of BMP 
percolation at low temperatures and RKKY-Zener cou- 
pling through activated carriers at high temperatures. 
Only detailed experimental work can establish the ori- 
gin of ferromagnetism in magnetically doped oxides — 
theory can only suggest interesting possibilities without 
ruling out alternative mechanisms such as extrinsic ferro- 
magnetism or direct superexchange. Among the various 
experimental evidence supporting the model of carrier- 
mediated ferromagnetism in mag netically doped oxides 
arc the observation of AHE llS], electric field induced 
modulation of magnetization |2l| . and optical magnetic 
circular dichroism 31|. It is clear that much more sys- 
tematic experimental data showing the magnetic proper- 
ties as a function of carrier density (and carrier proper- 
ties) will be needed before the issue of definitively estab- 
lishing the origin of ferromagnetism in doped magnetic 
oxides can be settled. 



V. CONCLUSION 

Before concluding, we emphasize that the new idea 
in this paper is that high-temperature RKKY ferromag- 
netism may be mediated in a semiconductor doped with 
magnetic impurities by thermally excited carriers in an 
otherwise-empty itinerant (conduction or valence) semi- 
conductor band, in contrast to the usual band-carrier me- 
diated RKKY ferromagnetism often discussed [H, 113] in 
the context of ferromagnetic (Ga,Mn)As where the va- 
lence band holes are thought to mediate the ferromag- 
netic RKKY coupling between the localized Mn mag- 
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netic moments. We believe that this RKKY coupHng 
mediated by thermally excited carriers may be playing 
a role in the high observed Tc in Co-doped Ti02 where 
the experimentally measured conductivity is always acti- 
vated in the ferromagnetic phase indicating the presence 
of substantial thermally activated free carriers. Obvi- 
ously, the necessary condition for such a thermally ac- 
tivated RKKY ferromagnetism is that electrical conduc- 
tion in the system must be activated (and therefore in- 
sulating) in nature in contrast to metallic temperature- 
independent conductivity. If the observed conductiv- 
ity in the system is temperature- independent, then the 
novel thermally-excited RKKY mechanism proposed by 
us simply does not apply. Our motivation for suggest- 
ing this rather unusual thermally-activated RKKY DMS 
ferromagnetism has been the reported existence of very 
high transition temperatures in Co-doped Ti02 which 
simply cannot be explained quantitatively by the bound 
magnetic polar qn p ercolation picture of Kaminski and 
Das Sarma lU IH| although it is certainly possible 
that some of the oxide ferromagnetism arises purely from 
the BMP mechanism as has recently been argued [s^l- 
The thermally activated RKKY ferromagnetic mecha- 
nism, while being necessary for high Tc, cannot be suf- 
ficient since at low temperatures the thermally activated 
carriers freeze out leading to an exponential suppression 
of the thermally activated RKKY ferromagnetism. We 
therefore propose that BMP ferromagnetism [i^, |43, 
mediated by strongly localized carriers take over at low 
temperatures, as already proposed by Coey and collabo- 
rators [lOl supplementing and complementing the high- 
temperature RKKY mechanism. The two mechanisms 
coexist in the same sample with the high values of Tc be- 
ing controlled by the thermally activated RKKY mech- 
anism and ferromagnetism persisting to T = due to 
the bound magnetic polaron percolation mechanism. We 
note that in a single sample the two mechanisms will 
give rise to a unique Tc depending on all the details 
of the system since the two mechanisms coexist whereas 



there could be considerable sample to sample Tc varia- 
tion due to the coexistence of the two mechanisms. In our 
picture the high-temperature thermally activated RKKY 
mechanism smoothly interpolates to the low-temperature 
BMP ferromagnetism with a single transition tempera- 
ture. This picture of the coexistence of two complemen- 
tary ferromagnetic mechanisms in oxides is essentially 
forced on us by our consideration of the possible tran- 
sition temperatures achievable within the BMP model, 
which are just too low to explain the observations in (at 
least) the Co-doped Ti02. 

Very recently, we have shown [7l| that the confluence 
of two competing FM mechanisms, namely the BMP per- 
colation in the impurity band at 'low' temperatures and 
the activated RKKY interaction in the conduction band 
at 'high' temperatures, could lead to an intriguing and 
highly non-trivial re-entrant FM transition in 0-DMS 
where lowering temperature at first leads to a non-FM 
phase which gives way to a lower temperature second 
FM phase. A direct observation of our predicted [7lj re- 
entrant FM would go a long way in validating the dual 
FM mechanism model introduced in this article. 

To summarize, we have analyzed different proposed 
models for carrier-mediated ferromagnetism in dilute 
magnetic oxides such as Ti02, ZnO, andSn02. Due to 
the insulating character of these compounds, a model 
based on the formation of bound magnetic polarons is 
proposed. However, the binding energy of the electrons 
on the oxygen vacancies that act as shallow donors is not 
large enough to keep the electrons bound up to the high 
temperature reported for the Tc (^ 700K). Therefore, 
we propose that, at sufficiently high temperatures, still 
below Tc, thermally excited carriers also mediate fer- 
romagnetism via an RKKY mechanism, complementing 
the bound polaron picture and allowing a considerable 
enhancement of Tc- 

We thank S. Ogale for valuable discussions. This work 
is supported by the NSF, US-ONR, and NRI-SWAN. 
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